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Критически проанализированы возможности и особенности существующих методов выявления и геометризации структур ано- 
мальных геохимических полей рудных месторождений. Установлено, что структура геохимического поля наиболее объективно 
выявляется в пространстве признаков. Предложены новые способы и методика комплексирования существующих методов ана- 
лиза структуры геохимических полей. 


Исследование зональности геохимических по- 
лей, как важнейшей составляющей комплексной 
зональности месторождений, изначально было од- 
ним из приоритетных направлений в прикладной 
геохимии. Являясь отражением последовательного 
отложения минеральных парагенезисов, геохими- 
ческая зональность была и остается ведущим фак- 
тором прогнозирования оруденения любого типа 
и, в особенности, гидротермального. 

После открытия в 70-х годах прошедшего века 
обобщенного ряда вертикальной геохимической 
зональности, закономерности последовательного 
отложения элементов стали использоваться во всех 
моделях зональности. Между тем, очень скоро вы- 
яснилось, что монотонное убывание с глубиной 
значений коэффициентов зональности часто соче- 
тается с очень сложным их распределением в плос- 
кости рудных тел, что делает проблематичным ис- 
пользование названных коэффициентов при дос- 


таточно редкой сети опробования. Особенно остро 
эта проблема встала для золоторудных месторожде- 
ний, где обычным является столбовое и ярусное 
распределение оруденения. Впоследствии в поня- 
тие аномалии, кроме ореолов привноса элементов, 
все чаще стали включать и области их выноса. Вы- 
яснилось также, что эти сопряженные в простран- 
стве положительные и отрицательные аномалии не 
являются чем-то изолированным, а представляют 
собой часть закономерно построенных конструк- 
ций со строгой иерархической соподчиненностью. 
Все это привело к появлению структурного метода 
исследования геохимических полей, при котором 
геохимические поля рассматриваются как структу- 
рированные системы, каждый блок которых харак- 
теризуется специфическими параметрами. В наи- 
более последовательном виде структурный метод 
изложен С. А. Григоровым [1]. Метод существенно 
облегчает разбраковку многочисленных аномалий, 
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выявляемых в процессе геохимических съемок, и 
позволяет выделить наиболее перспективные из 
них даже при идентичности количественных пока- 
зателей. Главной проблемой при этом является 
правильная интерпретация внутреннего строения 
геохимического поля, учитывая пространственное 
совмещение аномальных структур различных ие- 
рархических уровней. 

В настоящее время для расшифровки структуры 
геохимического поля используются, в основном, две 
группы методов: 1) выделение областей пространства 
со сходными геохимическими спектрами и 2) выяв- 
ление устойчивых ассоциаций элементов и анализ 
их пространственного размещения. Первая группа 
методов основана на кластеризации наблюдений, 
вторая осуществляет классификацию в признако- 
вом пространстве с последующей пространственной 
геометризацией выделенных ассоциаций. Наиболее 
широко используемые методы первого направления 
- метод многомерных полей, система "Геоскан" [2] , 
0-метод факторного анализа и его разновидности, 
кластер-анализ наблюдений, корреляция геохими- 
ческих спектров наблюдений, методики вычисле- 
ния показателей общей интенсивности рудного 
процесса (коэффициенты аномальности, комплекс- 
ный показатель дисперсии, ранговая дисперсия и 
др.), различные мультипликативные и аддитивные 
показатели зональности. Ко второй группе относят- 
ся К-метод факторного анализа, кластер-анализ пе- 
ременных, искусственные нейронные сети, дискри- 
минантный анализ, регрессионный анализ, другие 
методы классификации переменных. 

Определенные сведения о структуре геохими- 
ческого поля дают расчетные показатели, основан- 
ные на центробежно-центростремительной клас- 
сификации химических элементов Ю.Г. Щербако- 
ва [3]. Этим автором предложено рассчитывать ко- 
эффициенты относительной концентрации (ОК) 
родственных элементов, нормированные по не- 
дифференцированным хондритам. Максимальные 
значения этих показателей характерны для апи- 
кальных частей рудных тел, а в более мелком масш- 
табе, очевидно, - для флангов месторождений. 
В. И. Силаевым [4] предложены два показателя, ос- 
нованных на классификации Ю.Г. Щербакова, ко- 
торые мы предлагаем называть индексами центро- 
бежности (ИЦ1 и ИЦ2) и использовать их для ана- 
лиза структуры геохимического поля. Первый из 
них (ИЦ1) представляет собой отношение суммы 
минимально-центробежных элементов к сумме 
центростремительных элементов, второй (ИЦ2) - 
отношение суммы центробежных и дефицитно - 
центробежных элементов к сумме центростреми- 
тельных и минимально-центробежных. Содержа- 
ния элементов должны быть стандартизованы. 

Представляется полезным на конкретных при- 
мерах провести сравнительный анализ эффектив- 
ности перечисленных методик исследования внут- 
реннего строения геохимических полей различных 
иерархических уровней. Такой анализ проведен на- 


ми для разных типов золоторудных месторожде- 
ний, а также золотосодержащих скарново-магне- 
титовых и колчеданно-полиметаллических объек- 
тов. Ниже обсуждаются результаты по одному из 
золоторудных полей. 

Юхтинское рудное поле расположено в Цент- 
рально-Алданском золоторудном районе и вклю- 
чает в себя Самолазовское и Гарбузовское место- 
рождения и около десятка мелких проявлений и 
точек минерализации золота. Оруденение предс- 
тавлено субгоризонтальными залежами сульфид- 
но-флюорит-полевошпат-кварцевого состава и 
крутопадающими зонами прожилково-вкраплен- 
ной минерализации того же парагенезиса. Из суль- 
фидов, количество которых обычно не превышает 
5 %, развит, в основном, пирит, реже - галенит, 
халькопирит. Характерной особенностью орудене- 
ния рассматриваемого типа является его приуро- 
ченность к приконтактовой зоне штокообразных 
массивов сиенит-граносиенитов Лебединского 
комплекса (рис. 1). При этом рудные залежи лока- 
лизуются в скарнах, а прожилково-вкрапленная 
минерализация приурочена к зонам тектонических 
нарушений внутри интрузивных пород. 

В основу наших построений положены резуль- 
таты литогеохимической съемки масштаба 1:10000 
по вторичным ореолам рассеяния в пределах Юх- 
тинского массива и его обрамления (60 кв. км), вы- 
полненной ГУГПП "Алдангеология", и данные гео- 
химического опробования керна скважин, пробу- 
ренных на месторождении Самолазовском артелью 
старателей "Селигдар". 

Для всей площади работ во вторичном геохими- 
ческом поле факторным анализом выделяются сле- 
дующие устойчивые ассоциации главных элемен- 
тов-индикаторов: 1) РЬ, Ві, А§, Си; 2) Ва, Мп; 3) Аз, 
(А§, Аи), 4) Аи, Си, А§. Серебро и медь входят в сос- 
тав двух ассоциаций - золото-медной и свинцово- 
висмутовой. Характерно, что содержания меди и 
золота тесно коррелируются на участках с про- 
мышленным оруденением, на остальной площади 
медь входит только в состав первого фактора. 

На рис. 1, б показано пространственное разме- 
щение выделенных ассоциаций на изученной пло- 
щади. Они объединяются в субмеридиональную 
дугообразную полосу размерами порядка 10x3 км. 
В пределах этой структуры, по масштабу соответ- 
ствующей рудному полю, отчетливо выделяются 
3 зоны - центральная, сложенная перекрывающи- 
мися ассоциациями Аи, Си, А§ и РЬ, Ві, А§, Си, 
промежуточная, где преобладает ассоциация РЬ, 
Ві, Си, А§, и внешняя, с повышенными концентра- 
циями Ва и Мп. Промышленные месторождения 
обнаружены в центральной зоне, причем только 
там, где пространственно совмещаются две ассо- 
циации (факторы 1 и 4), в промежуточной и внеш- 
ней зонах известны только рудопроявления и точ- 
ки минерализации. 

Кластер-анализом наблюдений в рудном поле 
выявлено 4 устойчивых класса проб. Центральная 
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Рис. 1 . Модели геохимической зональности рудного поля по данным различных методов: А) схема геологического строения: 

1) доломиты; 2) мрамора; 3) песчаники; 4) граносиениты; 5) скарны; 6) основные разрывы; 7) золоторудные месторож- 
дения: 1 — Самолазовское, 2 — Г арбузовское; В) факторный анализ; участки развития ассоциаций: 1) РЬ, Ві, Ад, Си; 

2) Ва, Мп; 3) Аз, (Ад, Аи), 4) Аи, Си, Ад; С) кластер-анализ наблюдений; состав кластеров: 1) Ві, Си, РЬ, Аи, Ад; 2) фо- 
новые концентрации всех элементов; 3) Ва, Мп; 4) Аз, Ад, Си, Аи, РЬ; 0) индекс центробежное™ ИЦ2; Е) Г еоскан-мо- 
дель; Е) показатель ОКАд.Аи; С) ранговая дисперсия (по [5]); Н) энергия рудообразования 
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часть площади, с наиболее интенсивной минерали- 
зацией, отнесена к кластеру № 1 (Аи, Си, РЪ, Ві), к 
нему же отнесен слабо золотоносный участок на се- 
веро-западном фланге рудного поля (рис. 1, С). 
Кластер 4 (А§, А$) локализовал участок с предпола- 
гаемым оруденением куранахского типа. Кластером 
3, как и в факторном анализе, оконтурилась пло- 
щадь развития ассоциации Ва и Мп, остальные про- 
бы отнесены к фоновой совокупности (кластер 2). 
Таким образом, кластер-анализ наблюдений не смог 
дифференцировать участки совмещения минерали- 
заций различного типа в центральной части рудного 
поля. С учетом того, что интенсивность рудного 
процесса в результатах кластер-анализа наблюдений 
не отражается, очевидно, что результаты этого до- 
вольно трудоемкого метода во всех отношениях ус- 
тупают данным, полученным факторным анализом. 

Повышенные значения индекса центробежнос- 
ти ИЦ2=(А$+РЪ+Ві+Ва)/(Со+Ні+Си+А§+Аи) 
формируют кольцевую полузамкнутую структуру, 
окаймляющую участок максимального оруденения 
с севера, запада и юго-запада, фиксируя зону фрон- 
тального обогащения рудного поля (рис. 1, Л). 

На рис. 1, Е приведен результат исследования 
структуры геохимического поля Юхтинского мас- 
сива системой "Геоскан". Центральная часть рудно- 
го поля отнесена к кластеру РЬ 6 2 Ві 4 сАш^Сщ^А^і 7 , 
причем максимум развития этой ассоциации сов- 
падает с месторождением Гарбузовским. Самолазо- 
вское месторождение тоже находится внутри наз- 
ванного кластера, но интенсивность поля над ним 
достаточно низкая. Из сравнения с результатами 
факторного анализа видно, что на описываемой 
площади совмещены две геохимические ассоциа- 
ции. В результате кластером, где золото является 
ведущим элементом (класс № 7) "Геоскан" фикси- 
рует участки с непромышленной минерализацией 
на флангах известных месторождений. Тремя 
родственными классами проб медно-свинцово-се- 
ребряного состава с невысокой интенсивностью 
функции 8САК выделяются несколько локальных 
участков во фронтальной зоне рудного поля, на 
площадях с предполагаемым оруденением курана- 
хского типа. Таким образом, "Геоскан" уверенно 
выделил рудоносные площади, дифференцировав 
их по интенсивности рудного процесса и типу гео- 
химического спектра. В то же время, система не 
смогла идентифицировать участки перекрытия 
двух ассоциаций, медно-золотой и полиметалли- 
ческой, а именно эти участки и вмещают промыш- 
ленное оруденение. 

Аномалии коэффициента относительной кон- 
центрации ОК А§:Аи окаймляют месторождения 
Самолазовское и Гарбузовское с севера и с юга и, 
по всей видимости, маркируют рудоподводящие 
субмеридиональные структуры (рис. 1, Г). 

Максимальными значениями "показателей ин- 
тенсивности" (ранговая дисперсия и энергия рудооб- 
разования [5, 6]) уверенно выделяются оба промыш- 
ленных месторождения (рис. 1, С; 1, Я). Кроме того, 


на северном продолжении тектонической зоны, 
контролирующей размещение Гарбузовского место- 
рождения, локализуются еще 2 участка с аномально 
высокими значениями этих показателей. Промыш- 
ленных концентраций золота здесь не выявлено. 

В ранге месторождения по геохимическим дан- 
ным выделяется, как наиболее перспективный, 
участок, ограниченный ручьями Жильный, Карс- 
товый и рекой Большая Юхта (участок Жильный). 

Рис. 2, А демонстрирует структуру вторичного 
геохимического поля участка Жильного по данным 
факторного анализа. Здесь обособились 4 главных 
фактора, вклад которых в общую дисперсию сос- 
тавляет около 60 %. Они формируют очень контра- 
стную полузамкнутую концентрическую структуру, 
от центра к периферии которой последовательно 
сменяются ассоциации: (Мо, А§) - (Аи, Ак) - (РЬ, 
Ві, А§, 2п, Си) - (N1, Сг, Со). С юга и севера струк- 
тура окаймляется аномалиями Ва и Мп. Конфигу- 
рация аномальной структуры позволяет предпола- 
гать приуроченность ее к пересечению субмериди- 
ональных, субширотных и северо-восточных раз- 
рывов, из которых первые можно рассматривать 
как рудоподводящие. На участке локализуются 4 
аномальных структуры геохимического поля 
(АСГП), две из которых соответствуют месторож- 
дениям Самолазовскому и Гарбузовскому, две дру- 
гих фиксируют участки минерализации, связанные 
с северным продолжением субмеридиональных ру- 
доконтролирующих зон. Последние две структуры 
имеют неполный набор ассоциаций, что свиде- 
тельствует о недостаточной интенсивности рудно- 
го процесса. Промышленных концентраций золо- 
та, по данным буровых работ, здесь не выявлено. 

Системой "Геоскан" на участке Жильном выде- 
лено 7 аномальных кластеров с близкими геохими- 
ческими спектрами (рис. 2, В). Они концентриру- 
ются в нескольких локальных аномалиях, три из 
которых пространственно приурочены к Самола- 
зовскому и Гарбузовскому месторождениям, но 
фиксируют здесь участки развития полиметалли- 
ческих ассоциаций (РЪ, Ві, 2п, Си). Золото в них 
существенной роли не играет. Таким образом, оба 
промышленных месторождения системой "Геос- 
кан" выявляются, но структуру геохимического по- 
ля расшифровать по полученным данным без до- 
полнительной информации очень сложно. 

Кластер-анализом в пределах Жильного участ- 
ка выделено 5 классов проб (рис. 2, С). Золоторуд- 
ные месторождения уверенно фиксируются клас- 
тером № 2 с повышенными концентрациями РЪ, 
Ві, А§, 2п, Си, Т1, 8п, \Ѵ, V. Зоны непромышленной 
минерализации, включая участки на северном про- 
должении рудовмещающих зон, отнесены к клас- 
теру № 5, представленному теми же элементами, 
но с более низкими их концентрациями. В центре 
площади развита ассоциация Мо, А§ (кластер 
№ 4), а по периферии месторождений проявлен 
кластер № 3 (Со, №. Ті, Мп, Ій, Сг). Таким обра- 
зом, структура геохимического поля, выявленная 
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кластер-анализом наблюдений, близка к АГСП, 
полученной по результатам факторного анализа, за 
исключением интенсивности геохимических пре- 
образований, которая в результатах кластер-анали- 
за напрямую не отражается. 


Максимальные значения показателя энергии 
рудообразования характерны для Гарбузовского и 
Самолазовского месторождений, а также для двух 
участков на их северном продолжении. Последние, 
исходя из результатов кластер-анализа и факторно- 
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Рис. 2. Структуры геохимического поля участка Жильного. А)факторный анализ; участки развития ассоциаций: I) РЬ, Ві, Ад, Іп, Си; 
2) N1, Сг, Со; 3) Ва, Мп; 4) Мо, Ад; 5) Аи, Аз; 6) контуры месторождений: I — Самолазовского, II — Гарбузовского; 
В) Г еоскан-модель; С) кластер-анализ наблюдений; состав кластеров: I) фоновые концентрации всех элементов; 
2) РЬ 5 62 8/ 4і аз /Ад, -., 2п зд ТІ\ % Си <дз Аи уд 3) Со, ; 7 Л//і 05 77о 98 Мггд Сг а 84 ; 4) МОздАд 0 ^; 5) Сио ^Ащ^ в/ 0 ,б 7 77о,бо Адозі 2Пг,д РЬ 04 з; 
О) индексы центробежности: I) І4Ц1; 2) І4Ц2; Е) показатель ОКАд.Аи 
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го анализа, могут быть отнесены к зонам рассеян- 
ной минерализации. 

Повышенные значения индексов центробеж- 
ности в пределах участка Жильного приурочены к 
рудоконтролирующим субмеридиональным струк- 
турам, причем аномалии показателя ИЦ1 отчетливо 
окаймляют с двух сторон зону повышенных значе- 
ний показателя ИЦ2 (рис. 2, О). На месторождении 
Самолазовском это проявлено менее контрастно, 
чем на Гарбузовском, из-за пологого залегания руд- 
ного тела и значительной заболоченности участка. 

Высокие значения коэффициентов относитель- 
ной концентрации пары элементов А§:Аи приуроче- 
ны к тем же рудоконтролирующим структурам и ха- 
рактерны для флангов месторождений (рис. 2, Е). На- 
иболее золотоносные участки характеризуются пони- 
женными значениями ОК и ИЦ. По этому признаку 
северное продолжение Гарбузовского месторождения 
может быть классифицировано как зона рассеянной 
минерализации, что подтверждает выводы, получен- 
ные кластер-анализом и факторным анализом. 

Зональность размещения геохимических ассо- 
циаций в плоскости рудных залежей рассмотрена 
на примере разведочных пересечений через Само- 
лазовское месторождение (рис. 3). По результатам 
факторного анализа в пределах пологозалегающей 
залежи выявлено 4 класса элементов. Наиболее зо- 
лотоносная часть залежи характеризуется аномаль- 
ным накоплением Аи, А§, Аз, Т1, РЬ, Си, Мо, №, V, 


8Ь, \У. При этом Аз, N1, Т1, 8Ъ тесно коррелируются 
с золотом и объединяются с ним в одну ассоциа- 
цию. Участки ее распространения в целом совпада- 
ют с контурами промышленных блоков. Ассоциа- 
ция РЬ, Ві, А§, 2п, Си распространена более широ- 
ко, но преимущественно в восточной части место- 
рождения. В западном направлении она сменяется 
на ассоциацию Ва, Мп, Сг, которая сопровождает- 
ся умеренной и слабой золотоносностью. Харак- 
терно, что наиболее золотоносные участки сопро- 
вождаются локальными аномалиями титана в ле- 
жачем и висячем боках рудных тел. Ті коррелирует- 
ся с Мо и, в целом, эта ассоциация распространена 
на участке в виде полосы северо-восточного прос- 
тирания. Возможно, она фиксирует рудоконтроли- 
рующую структуру ранга месторождения. 

Системой Теоскан" в разрезе через рудную за- 
лежь выявлено 7 аномальных классов проб со сход- 
ными геохимическими спектрами, в каждом из ко- 
торых главную роль играет РЬ. Близость состава ге- 
охимических спектров, обусловленная простран- 
ственным совмещением различных геохимических 
ассоциаций, не способствует пониманию механиз- 
ма формирования геохимической зональности. 
Промышленная часть залежи выделяется, по-су- 
ществу, лишь повышенным вкладом Аи в геохими- 
ческий спектр (рис. 3, В). Структуру же геохими- 
ческого поля, без априорной информации, полу- 
ченной другими методами, расшифровать доста- 



Рис. 3 . Структура геохимического поля ранга рудного тела. Разрез через месторождение Самолазовское: А) факторный анализ: 

1) граносиениты ; 2) дезинтегрированные скарны с рудной минерализацией; участки развития геохимических ассоциа- 
ций: 3) Аи, Аз, N1, ТІ, 5Ь; 4) РЬ, Ві, Ад, 2п, Си; 5) Ва, Мп, Сг; 6) Ті, Мо; 7) разведочные скважины; В) кластер-анализ наб- 
людений; состав кластеров: 1) Віоп РЫ.^ СГо ,48 Со$, & 2) Г/о,з 8 Ве$ 28 Мп^я Ваг.-.-.; 3) Цо? А/г. с, Г (, 51 Ад-. А] ЗЬ,^Аи.^Аз-.,;.і РЬ,^. Сц 44 МОі&; 
4) Г /,,05 &э 0 і іб,' С) области повышенных значений индексов центробежности: 1) ИЦ1=(Р+ТІ+Ѵ)/ (Мп+Си+Ад+Іп+ЫІ+Со+Аи; 

2) ИЦ2=(РЬ+Аз+5Ь+Ва+ВІ+Мо)/(Р+ТІ+Ѵ+Мп+Сг+Си+Іп+Ад+ЫІ+Со+Аи); 3) рудное тело; й) участки повышенных значе- 
ний показателя ОК Ад.Аи; Е) Г еоскан-модель 
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точно сложно. Следует заметить, что в этих услови- 
ях обычный кластер-анализ наблюдений по нор- 
мированным данным дает более наглядную карти- 
ну зональности (рис. 3, С). 

Весьма показательно распределение в плоскости 
разреза значений индексов центробежности (рис. 3, 
/)). Концентрический характер распределения зна- 
чений ИЦ1 и ИЦ2 вокруг рудного тела дополняет 
картину зонального строения АСГП, выявленную 
вышеописанными методами, и позволяет говорить 
о принципиальной возможности создания типовой 
геохимической модели подобных объектов. 

Распределение значений О К А§:Аи в плоскости 
разреза свидетельствует о приуроченности их мак- 
симумов к слабозолотоносной периферии рудной 
залежи (рис. 3, Е). Таким образом, данный коэф- 
фициент может быть использован как для просле- 
живания рудовмещающих структур, так и для рас- 
шифровки структуры геохимического поля. 

Повышенные значения показателя энергии ру- 
дообразования однозначно коррелируютя с содер- 
жаниями золота. Никакой новой информации этот 
коэффициент не дает, поэтому его использование в 
масштабах первичных ореолов отдельных рудных 
тел, видимо, излишне. 

Итак, характерными геохимическими ассоциа- 
циями самолазовского типа оруденения являются 
1) Аи, Аз, №, Т1, 5Ъ; 2) РЬ, Ві, А§, Си, 2п; 3) Ва, Мп. 
Состав их несколько видоизменяется на различных 
иерархических уровнях вследствие частичного пе- 
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рекрытия, но в целом выдерживается достаточно 
хорошо. Характерно, что АСГП различных уров- 
ней, от рудного поля до рудных тел, являются фи- 
гурами подобия, особенностью их внутреннего 
строения является приуроченность Аи и ассоци- 
ирующих с ним элементов к центрам аномальных 
геохимических структур. 

Внутреннее строение АСГП на всех иерархичес- 
ких уровнях наиболее полно и объективно выявля- 
ется с использованием методов кластеризации пе- 
ременных, лучшим из которых является фактор- 
ный анализ. Дополнительная информацию дают 
кластер-анализ наблюдений, система Теоскан", 
показатели интенсивности рудного процесса, при- 
чем роль этих вспомогательных методов снижается 
с увеличением детальности исследований (от руд- 
ного поля к рудному телу). Напротив, показатели, 
основанные на центробежно-центростремитель- 
ных свойствах элементов (ОК, ИЦ1, ИЦ2), наибо- 
лее интересные результаты дают в масштабе руд- 
ных тел. Сходные данные получены и для других 
типов изученных месторождений. Комплексирова- 
ние перечисленных методов, с учетом обоснован- 
ной выше приоритетности, позволяет создавать 
обобщенные геохимические модели рудных объек- 
тов. Подобные морфологические модели могут 
быть описаны количественно и типизированы 
применительно к различным группам рудных мес- 
торождений и различным уровням их эрозионного 
среза. Создание таких типовых моделей - насущ- 
ная проблема сегодняшнего дня. 
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